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ABSTRACT 

Using a new Chandra ^C/S'-S3 observation of the normal, isolated, moderate-Lx, lenticular galaxy 
NGC 1332 we resolve the emission into ~75 point-sources, and a significant diffuse component. We 
present a detailed analysis of the spectral properties of the diffuse emission, constraining both the 
temperature profile and the metal abundances in the hot gas. The characteristics of the point source 
population and the spatial properties of the diffuse emission are discussed in two companion papers. 
The diffuse component comprises hot gas, with an '^isothermal temperature profile (^0.5 keV), and 
emission from unresolved point-sources. In contrast with the cool cores of many groups and clusters, 
we find a small central temperature peak. We obtain emission-weighted abundance contraints within 
20 kpc for several key elements: Fe, O, Ne, Mg and Si. The measured iron abundance (ZFe=l-l in 
solar units; >0.53 at 99% confidence) strongly excludes the very sub-solar values often historically 
reported for early-type galaxies. This continues, in a lower-Lx system, a trend in recent observations 
of bright galaxies and groups. The abundance ratios, with respect to Fe, of the other elements were 
also found to be ~solar, with the exception of Zo/Zpe which was significantly lower (< 0.4), as seen in 
several bright galaxies, groups and clusters. Such a low O abundance is not predicted by simple models 
of ISM enrichment by Type la and Type II supernovae, and may indicate a significant contribution 
from primordial hypernovae. Revisiting Chandra observations of the moderate-Lx, isolated elliptical 
NGC 720, we obtain similar abundance constraints (ZFe=0.7lto;2iJ 90'^ confidence; Zo/ZFe=0.23 ± 
0.21). Adopting standard SNIa and SNII metal yield models, our abundance ratio constraints imply 
73±5% and 85±6% of the Fe enrichment in NGC 1332 and NGC 720, respectively, arises from SNIa. 
Although these results are sensitive to the considerable systematic uncertainty in the SNe yields, they 
are in good agreement with observations of more massive systems. These two cases of moderate-Lx 
early-type galaxies reveal a consistent pattern of metal enrichment from cluster scales to moderate 
Lx/Lb galaxies. 

Subject headings: Xrays: galaxies — galaxies: elliptical and lenticular, cD — galaxies: abundances — 
galaxies: individual(NGC 1332) — galaxies: halos — galaxies: ISM 



1. INTRODUCTION 

During the formation and evolution of early-type 
galaxies, the primordial gas is enriched by supernovae 
ejecta and stellar mass-loss. Observations of the metal 
content of the hot, X-ray emitting gas within such galax- 
ies therefore offer a powerful means to inves tigate their 
histo ry (e.g. Loowonotoin & Mathowc - 1901 )-. — In early 



the giant elliptical galaxies (Renzini 1997). This discrep- 
ancy is problematical, as it is difhcult to envisage how 
individual early- type galaxies could have, in many cases, 
lower iron abundances than typical clusters. Indeed, gas 
enrichment models for the centres of individ ual galaxies 



tend t o predict super-solar values of Zf^ fe.g. Ciotti et al 
Loewenstein & Mathews 1991; Brighenti fc MatE 



1991 



ASQA and Rosai observations of earljf-type galaxies vcr}^ 
sub-solar metal abundances (dominated by Fe, which has 
the strongest diagnostic lines in the soft X-rav band) were 



^ 1999| ) . This obvious inconsistency led [Arimoto et al 



(|1997) to question the validity of the X-ray plasma codes 



j:enerally reported (e.g. Davis &: White _ 199 




which were being fitted to the data, especially given un- 
certainties in the physics of the Fe L-shell tran sitions. 



An a lt ernat ive explanation was provided by 



Buote 



Fabian (1998), who demonstrated that fitting a single 



contrast, the mean stellar ir on abundance in elli ptical 
galaxies tends to be ~solar ( Arimoto et al. |1997 ), im- 
plying that the ISM could not have been substantially 
enriched by SNIa and mass-loss from the stellar popu- 
lation. However, the mean iron abundances determined 
by X-ray observations of clusters tend to be ~0.3-0.5 
times solar, requiring significant enrichment of the pri- 
mordial gas, most probably from the stellar population of 
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temperature model to intrinsically non-isothermal data 
can give rise to an "Fe bias" , in which the iron abundance 
is systematically under-estimated, and instead obtained 
~solar abundances by fitting multi-temperature models 
to a v arietv of bright elliptical gala xies (see also Buote 
2000b and e.g. Buote et al, 2003b, who found compa- 
rable results with different plasma models). Although it 
had long been recognized that a hard spectral component 
(to account for unresolved point source emission) was re- 
uired to fit the ^6' C^ data of many early-type galaxies 
[Matsushita et al. 1994), the spectral shape of this com- 



ponent is much harder (kT> 5 keV) than the hot gas 
emission (for which kT'-~^0. 5-1.0 keV), so that it cannot 
move to mitigate, even in part, the Fe bias. Nonetheless, 
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Buote fc Fabian ( 1998 ) still found that, in the lower-Lx 
systems, in which the data were consistent with a single 
hot gas component plus unresolved sources, the abun- 
dances were consistent with Zpc— 1- However, the lower 
signal-to-noise resulted in very poor constraints so that 
very sub-solar Zpc values could often not be entirely ruled 
out . A similar result for a low S/N system was reported 
by |Kim et al.| dl99d ) for the low Lx/Lb SO NGC4382 
observed with ASCA. 

More recent Chandra and XMM observations of bright 
ellipticals and the centres of bright groups have tended to 
confirm ne ar solar (or even super-solar) abundanc es (e.g 
NGC 5044:|Buote ct al | | 2003b| ,[ramur a et al] |2003|: M87 



pastaldello fc Molcnd i "p Op 
toj4636: ku ct a[\\2m: 



NGC 13 99nBuotL] |2002 
'MKW4: IO 'Sullivan eOf 



003; also c.f. the low Lx/Lr radio galaxy NGC 1316 
im fc Fabbianc 2003). However, for some systems, es- 
pecially (but not exclusively) the lowest-Lx/LB galax- 
ies, authors are still tending to fi nd very sub-solar 



metal abundances 

NGC 4697: paralin et al 

NGC 4494 and NGC 5322 



e.g. NGC 12 91 



2001 
for 



2002, 



bright radio galaxy NGC 6251: ^ambruna ct al. 2004) 



nl 
t 



Irwin et al, 
especially JNGC 3585 
which Zpo^JO.l was re- 



ported by 0' Sullivan fc Ponmaij 2004|: also c.f. the X-r ay 



The lack of a consistent picture of enrichment from 
galaxy to cluster scales is a major problem in our un- 
derstanding of galaxy evolution. It is therefore of crit- 
ical importance to determine accurate abundances in 
galaxies with a wide range of Lx/Lb, and especially 
to determine the ubiquity of the very sub-solar abun- 
dances in 1ow-Lx/Lb systems. In this paper, we present 
a Chandra study of the metal abundances in the hot 
gas of the "normal", relatively isolated, moderate-Lx 
l enticular galaxy NGC 133 2 . In two compani on papers 
( [Humphrey fc Buot^ |2004| ; [Buote et al| ^004 hereafter 
Paper I and Paper III), we present a study of the X-ray 
point-source population and the gravitating matter dis- 
tribution, which have been able to take advantage of the 
excellent spatial resolution of Chandra. 

NGC 1332 has a much lower Lx (and Lx/Lb) than 
most of the gas-rich giant elliptical galaxies for which 
^solar abundances have been found. Furthermore, its 
Lx/Lb is sufficiently high to ensure substantial emission 
from the hot gas. For many systems with much lower 
Lx/Lb, such analysis is hampered by a lack of pho- 
tons (which may introduce abu ndance degeneracies o r 
extremely poor constraints; c.f. Buote fc Fabian 1998). 
Significant diffuse gas emission in NGC 1332 was inferred 
from previous X-ray observations, making it an ideal sub- 
ject for such a study. Us ing Rosat PSPC observations. 



Buote fc Canizares (1996) found evidence of a flattened 
dark matter halo in NGU1332. Although the ^S'C^ data 
enabled the unambiguous identification of a hard spec- 
tral component (supposed to arise from unresolved point 
sources), the signal-to- noise was insufficient to enable 
good c onstraints on the met al abundances to be deter- 
mined ( Buote fc Fabian] 1998[ ) (rcscaling to the new solar 
abundances, these authors found Zfo> 0.3). Although 
the best-fit value (Zpe— 1-2) was clearly consistent with 
^solar abundances, highly subsolar values could not be 
excluded. 

Throughout this paper, we adopt th e new standard 
abundances of Grevesse fc Sauval ( 1998| ). We note that 
there is some complication in comparing to the literature 



as many authors still u s e the outdated solar abundances 
of Anders fc Grevesse (1989). Although we quote all 
abundances with respect to the new standard, the most 
significant discrepancy between fits using the two stan- 
dards is a difference in the inferred value of Zpo; for com- 
parison with our work, Zpo determined using the older 
abundances should be multiplied by ^1.5. All quoted 
confidence regions are 90%, unless otherwise stated. 

2. BACKGROUND AND FLARING 

The region of sky containing NGC 1332 was observed 
with the ACIS instrument aboard Chandra between 2002 
September 19 10:39 and September 20 02:59 UTC, the 
galaxy being centred on the S3 chip, for a nominal 
~60 ks exposure. An additional observation made be- 
tween September 18 02:56 and 08:27 UTC, for a nom- 
inal 20 ks exposure, was heavily contaminated by flar- 
ing of the background and so we present results only 
for the better-quality dataset. We reduced and pro- 
cessed the data using the CI AG data-analysis suite of 
tools, version 3.0.2, and the Chandra Caldh version 2.26, 
and we performed spectral-fltting with XSPEC version 
11.3.0. For a discussion of the data-reduction procedures 
and a detailed analysis of the point-source population in 
NGC 1332, we refer the interested reader to Paper I. 

One of the challenges in performing analysis of rela- 
tively low surface-brightness diffuse emission is the deter- 
mination of a reliable background estimate. Since the dif- 
fuse emission covers much of the S3 chip, it was not pos- 
sible to determine a "local" background estimate from an 
appropriate region of the detector. We therefore adopted 
the appropriate "blank fleld" event flies provided in the 
Caldb. There is known to be variation between different 
observations and different pointings in both the high and 
low-energy regimes of the background. This is compli- 
cated further by the "flaring" which is often observed in 



the background (e.g. Markevitch 2002). We therefore ac- 
cumulated a background lightcurve (from 5.0-10.0 keV) 
from source-free regions of the active chips and examined 
it by eye for evidence of flaring. The lightcurve shows ev- 
idence of mild flaring activity. The strongest ffares had 
amplitudes reaching ~50%, but lasted for only a short 
fraction of the entire observation. In contrast, low-level 
flaring persisted through much of the observation, with 
an associated amplitude of ^20%. 

Completely excising the low-level flaring from the data 
resulted in a total exposure of only ~28 ks. We found 
excellent agreement in the high-energy count-rates of the 
background templates and the flare-cleaned observation, 
indicating complete removal of the flare. However, such a 
short exposure unacceptably degraded the signal-to-noise 
level of the data. As a compromise, we excised only the 
periods of strongest flaring, which resulted in a more ac- 
ceptable ~45 ks exposure. We chose not to include all of 
the data since the signiflcantly higher background during 
the remaining time periods would more than offset the 
signal-to- noise gain from the modest increase in exposure 
that would result. However, we compared spectra accu- 
mulated with each level of filtering from the same region 
of the S3 chip containing low surface-brightness emission 
from NGC 1332. The less- filtered data clearly showed 
significant residuals in excess of the hot gas emission at 
> 2 keV which were not seen in the more heavily-filtered 
data. In order to take account of this contamination, we 
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Fig. 1. — The spectra of the flaring, shown along with the best-fit models. Left panel: the flare spectrum obtained from the S3 chip. 
Right panel, the flare spectrum obtained from the SI chip. Note that, since the extraction regions were not identical, we do not expect the 
spectral normalizations to be identical. 



chose to include an additional component in our spectral- 
modelling. We did not attempt to renormalize the back- 
ground template to match the high-energy residuals as 
this will also renormalize the soft X-ray background and 
instrumental features (which can be problematic if the 
flaring is sufficiently strong, as in Sect |). 

To enable us to estimate the spectralform of the flar- 
ing, it was necessary to disentangle it from the diffuse 
galactic emission and the X-ray background. To achieve 
this, we adopted two complementary techniques. Firstly, 
we exploited the strong variability of the flaring compo- 
nent, in contrast to (naturally) constant galactic emis- 
sion. We accumulated a spectrum at the peak of flaring 
and one during a period when flaring was absent, from 
the same annular region from 5-160" (excluding all pho- 
tons within a region 6 times the radius of the la de- 
tection region around each resolved point-source, since 
individual sources may be variable). The difference be- 
tween these two spectra must, necessarily, be the spec- 
trum of the flare. The two periods were chosen from the 
lightcurve by eye (we note that intensity-filters would 
not have been appropriate since the Poisson scatter, for 
realistic lightcurve binsizes, was not much less than the 
flaring amplitude) . Appropriate count-weighted spectral 
response matrices were generated for the peak flaring 
data with the CIAO tasks mkrmf and mkvirarf, and 
the flare-free spectrum was adopted as the background. 
The "flare" data were regrouped to ensure at least 20 
source-plus-background photons in each spectral bin, and 
a signal-to-noise ratio of at least 3. This procedure as- 
sumes that the time-averaged spectrum of the flare is 
similar to the peak flare spectrum. A comparison of flare 
lightcurves accumulated in different energy bands (0.5- 
3.0 keV and 3.0-10.0 keV), obtained by subtracting the 
"quiescent" level (estimated from non-flaring stretches 
of data) from the background lightcurves, did not re- 
veal any significant changes in spectral hardness during 
evolution of the flare. This would seem to support our 
assumption. 

To parameterize the spectral shape of the flare we 
tested various spectral models against the data (which 
were not required to be physical). We found that the 
data could be very well-fitted (% /dof = 41.5/46) with 



a broken-power law model, plus a broad gaussian term 
at '^2.6 keV to soften the change in the spectral slope. 
In order to determine whether the flare spectrum varied 
with position on the S3 chip, we extracted similar spectra 
from the set of 5 annuli adopted in Sect 0. Excluding the 
two innermost annuli (for which the data were too poor 
for comparison), we found that the parameterized model 
gave an excellent flt to the data within each annulus, if 
the overall normalization was allowed to fit freely. The 
normalization scaled approximately with the extraction 
area, as expected for flaring which was uniform across 
the S3 chip. 

Our alternative means of estimating the flare spectrum 
was to accumulate the spectrum from a source- free region 
of the SI chip (since this was sufficiently offset from the 
galaxy centroid not to be strongly contaminated by emis- 
sion from the hot gas) and account for the X-ray back- 
ground by adopting the appropriate background tem- 
plate. We accumulated a spectrum from a region centred 
on the SI chip, with a 160" radius. We excluded data in 
the vicinity of a single point-source found by the detec- 
tion algorithm on this chip. We extracted a background 
spectrum from the standard SI "blank-fields" events file. 
The source spectrum was grouped to ensure there were 
at least 20 photons in each data bin and the signal-to- 
noise ratio exceeded 3. We were able to parameterize the 
SI flare spectrum with a single broken power law model. 
There was excellent agreement between flare spectra ac- 
cumulated via both methods above ~2 keV (Fig^). Be- 
tween ~l-2 keV, the flare accumulated from the si chip 
showed a significantly higher count-rate than that ob- 
tained from the S3 chip. A number of systematic effects 
may account for this discrepancy. Firstly, the SI chip 
may not be entirely free of contamination by the diffuse 
emission from the galaxy. Secondly, the degradation of 
the PSF on the SI due to the large offset means that 
point-sources may not have been detected on this chip, 
contaminating the data. Thirdly, the actual flare spec- 
trum may be different between the two chips. Finally, 
there may be subtle differences between the background 
template and the actual background spectrum for our 
pointing. In the subsequent sections, we adopt the flar- 
ing model taken from the S3 chip, since, on balance, it 
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is less likely to be subject to biases than the data from 
the SI chip. However, we found very little diff erence if 
the alternative flare model was used (Sect [4.4|) . Given 
the very different assumptions involved in our two esti- 
mates of the flare spectrum, this gave us confidence in 
our ability to model the flaring. 

3. SPECTRAL ANALYSIS 



Fo ■ analysis of the diffuse galactic emission of & Lockman 



excellent fit to the composite spectrum of the resolved 
sources (Paper I). To model the flaring, we used the pa- 
rameterization adopted in Sect g for the flaring on the S3 
chip, albeit with the normalization free. The components 
representing emission from the galaxy were additionally 
modified by photoelectric absorption due to cold material 
in the line-of-sight. The absorbing column was f ixed at 
the Galactic valu e for the appropriate pointing (Dickey 



NGG l'6'62, we extracted spectra from a smgle, large (3') 
aper ture and a series of contiguous, concentric annuli, 



sections of V 



L99C ) , and we adopted the ab sorption cross- 



orrison & McGammon (1 983, but the results 



all centred on the centroid of the X-ray emission. We 
determined the galaxy centroid self-consistently by it- 
eratively computing the centroid in a 1' circle, which 
was re-centred at this new position, until no more re- 
finements were required. Our initial guess for the cen- 
troid of X-ray emission was the 2MASS position given in 
NED. Our resulting position was in excellent agreement 
with that of the diffuse source identified at the galaxy 
centre by the point-source detection algorithm (Paper I) 
and was within 1" of the infra-red centroid. In order to 
minimize contamination from point sources we excluded 
regions around each detected point source (Paper I) of 
radius 6 times the semimajor axis of the l-cr encircled 
energy ellipse. The size of the large aperture was chosen 
to maximize signal-to-noise, whilst the individual annuli 
were chosen so as to contain in each approximately ^^1500 
background-subtracted photons. This gave annuli with 
outer radii 0.07', 0.55', 1.4', 2.3' and 3.1'. Appropri- 
ate count-weighted response matrices were generated for 
each spectrum, using the standard CIAO tools mkrmf 
and mkwarf. The spectra within each annulus were re- 
binned to ensure at least 20 photons in each spectral bin 
and to ensure a signal-to-noise ratio of at least 3. 

We restricted our fitting to the energy-band 0.5- 
7.0 keV, on account of calibration uncertainties at lower 
energies and the steeply-rising background at high ener- 
gies. At low energies, the background template spectrum 
becomes less reliable since there is known to be strong 
dependence in the soft X-ray background emission with 
pointing. Wc therefore took care to ensure that the data 
were truncated at low energies where the background and 
source began to converge. For the spectra of the two out- 
ermost annuli, this involved truncating the data below 
~0.7 keV; for the innermost annuli data were excluded 
below 0.5 keV. Gare should be exercised when compar- 
ing results obtained for different energy bandwidths, as 
narrowing the bandwidth can introduce systematic un- 
certainties in the fitting. However, since the abundances 
were subsequently tied together in all annuli, we found 
that they were relatively unaffected by this procedure 
(Sect 4.5). We discuss the likely impact of uncertainties 
in the background in Sect 42 . 



were unaffected if we adopted those of Batucihska-Church 



fc McGammon 1992) 



3.1. The large- aperture spectrum 

We modelled the spectrum with a hot gas component, a 
term to account for unresolved point-source emission and 
a component to accommodate the mild flarin g. To model 



t he ho t gas, we adopted the APEC model of Smith et al 



( ^OOlD 



since it was the most up-to-date plasma code 
which is widely available (the impa ct of using a MEKAL 
model instead is addressed in Sect O). To account for 
unresolved point-sources, we included a bremsstrahlung 
model with kT = 7.3 keV, which was found to give an 



We allowed Zpc to be free and initially tied the abun- 
dances of all remaining elements to Fe in their appropri- 
ate solar ratios. This model gave a good fit to the data 
(xVdof=121.04/114), with Zfo> 1.5. However, there 
were obvious sytematic discrepancies between the model 
and the data at low energies (^0.7 keV). We experi- 
mented with adding an additional hot gas component at 
a different temperature and, alternatively, systematically 
freeing individual metal abundances. In practice, the 
low-energy residuals were only significantly reduced by 
freeing O, for which the improvement in x^ was > 99.9% 
significant, on the basis of an F-test. We obtained best- 
fit abundances of ZFc=0.77ta32 and Zo/Zfc=0.23J:oJ3, 
for a gas temperature of 0.57±0.02 keV. The fit was ex- 
cellent (x^/dof=106/113). We experimented with sys- 
tematically freeing the remaining abundances. Although 
there was no significant improvement in the fit statistic 
following this procedure, we were able to estimate rea- 
sonable 90% confidence limits for the abundances of Ne, 
Mg, and Si (Table mf". The spectrum and best-fit model 
are shown in Fig g. The lines of Ne, Mg and Si are 
clearly visible in the spectrum, and so it is unsurprising 
that we have been able to constrain their abundances. 
In order to minimize the number of variable parameters, 
we kept fixed the shapes of the flare and point-source 
model components. To determine the impact of this as- 
sumption on our error-bars, we additionally fitted the 
3' aperture spectrum simultaneously with the flare and 
composite point-source spectra. Both the best-flt values 
and the error-bars were minimally affected by this proce- 
dure. We further investigate the impact of varying these 
components in Sect H. 

Although the low-energy residuals were substantially 
reduced by freeing Zq, there was still some evidence of 
a feature at these energies (although, due to the good 
flt at higher energies, it was not reflected transparently 
in a poor overall y^ value). We found that adding an 
additional hot gas component (which we migh t ex pect 
given the temperature gradient observed in Sect 3.2) led 
to a statistically significant improvement in fit (at 99% 
significance, on the basis of an F-test). The quality of 
the fit was excellent (x^/dof=97/lll). The tempera- 
tures of the two hot-gas components were found to be 
0.33lo°4 and 0.72to:?5 keV (in good agreement with the 
range of tem peratures seen over the temperature gradi- 
ent; Sect ^.2| ). The best-fitting abundances are shown in 
Table |^, and clearly indicate ~solar Fe abundance, and 

■^ Note: since the abundances ratios of Ne, Mg and Si to Fe 
were all consistent with unity, we determined the confidence region 
for each one of these elements assuming the abundance ratios for 
the others were exactly unity. This step was necessary to obtain 
interesting constraints 
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Fig. 2. — Left panel: spectrum of the large aperture. Right panel: spectrum of the innermost annulus. Also shown are the best-fit 
spectral-models (solid line; we use the single-temperature fit for the large aperture) and the individual contributions due to the hot gas 
(dotted line), the "contamination" due to projection effects (dashed line, right panel only), the unresolved point-sources (dash-dot line) 
and flaring (dash-dot-dot). The significance of the flaring component rises in the outermost annuli. Note that the feature at ~0.72 keV in 
the left panel is too narrow to be an emission line; it is a statistical fluctuation at < 3cr in a single data-bin. Since these spectra are not 
independent, it is visible in both panels, but it is not seen in other annuli. 



TABLE 1 

Emission-weighted abundances within the single 3' (20 kpc) 

aperture 



Model Zpo 


Zo/Zpo 


Znc/Zfc 


ZMg/Zpo 


Zsi/Zpc 


0.5-7.0 keV 


IT 0.77t,V.°32 
2T 1.09(>0.60) 


^■^■^-0.23 
0.02(<0.20) 


-, r>-|+U.bb 

t.^i_0.41 


1.08±0.31 


1.08±0.59 


0.4-7.0 keV 


IT i.iot^i? 

2T >1.03 


0.34t;j;-^ 
0.15 ±0.14 


-1 r)-|+U.b3 

J^-^J^-0.41 


1.04±0.31 


1.02±0.57 



Note. — Abundance measurements and 90% confldonce intervals de- 
termined from fitting the large-aperture spectrum. We present both single 
temperature (IT) and two-temperature (2T) results for the bandpass 0.5— 
7.0 keV and 0.4-7.0 keV. 



significantly sub-solar O. We experimented with allow- 
ing the remaining element abundances to vary, as for 
the single-temperature case. However, we found that 
none of the abundances (excepting Zpe and Zq) could 
be constrained. Since adopting an over-simplistic spec- 
tral model can give rise to spurious abundances, we ad- 
di tionally tested the " power law DEM" (PLDEM) model 
of Buote et al. (2003a), in which the differential emission 
measure is a power law function of the temperature, and 
the power law exponent and range of gas temperatures 
are determined as fit parameters. The inferred abudances 
(Zfc>1-1) and (Zo/Zfc< 0.22) were fully consistent with 
our two-temperature re sults. This agreern ent. which was 
also seen in NGC 5044 ( Buotc ct al. 2003b| ), gives us con- 
fidence in the validity of the abundances derived from 
our two-temperature parameterization of the spectrum. 
This is further supported by the good agreement with 
the spatially-resolved analysis, below. 

The ability to obtain reliable abundance measurements 
is critically dependent on our continuum emission con- 
straints. Especially when there are multiple spectral 
components (i.e. hot gas, unresolved point sources and 
a flare component), it is essential to constrain the con- 
tinuum at both ends of the spectrum. Below ^0.5 keV, 



the spectrum is relatively free from strong emission lines, 
and so the data are extremely useful to this end. Un- 
fortunately, the calibration of the ACIS chips is more 
uncertain in this regime, in part due to the build-up of 
carbanaceous deposits on the optical filter. Nevertheless, 
since we were using the most up-to-date correction for 
this quantum-efficiency degradation, we experimented 
with extending the bandpass down to 0.4 keV. Using ei- 
ther 1 or 2 temperature best-fit model for the 0.5-7.0 keV 
bandpass, we noticed a slight, systematic over-estimation 
of the count-rate below 0.5 keV (the discrepancy being 
more pronounced for the 2 temperature case). Fitting 
the data in the extended band, we found systematically 
higher values of Zpe and Zo/Zpe for each model (Table |^), 
but the remaining abundance ratios were relatively un- 
affected. Although there are calibration uncertainties in 
this regime, we observed very similar behaviour using 
an AOl observation of NGC 720 (Sect |), for which the 
ACIS contamination was far less severe, which gave us 
confidence in these results. 

3.2. Spatially-resolved spectroscopy 

Observations of giant galaxies and groups often iden- 
tify strong temperature gradients. In order to investigate 
any such effects, we performed spatially-resolved spec- 
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troscopy in a series of concentric annuli, centred on the 
peak of the X-ray emission. If there is a significant tem- 
perature gradient in the data, on account of the quasi- 
elHpsoidal distribution of the hot gas in the galaxy halo, 
the projection of the spectrum onto the sky does not ac- 
curately reflect the physical conditions at any point in the 
gas. In order to account for this effect, we employed the 
"o nion-p c cling" spectral deprojection algorithm outlined 
in Buote (2000a). To determine the extent to which such 
projection effects may alter our results, we additionally 
fitted traditional models directly to the spectra projected 
onto the sky. Henceforth we distinguish between pro- 
jected and deprojected results by the soubriquets "2D" 
and "3D", respectively. 

The quality of our data did not enable us to deter- 
mine the metal abundance independently in each annu- 
lus, preventing our investigation of any abundance gra- 
dients. For the 2D fitting, to maximize our constraints 
on the data, we fitted the projected spectra simultane- 
ously and tied equivalent metal abundances between each 
annulus. Since the magnitude of flaring did not vary 
strongly with position on the S3 chip (Sect ||), we were 
able to constrain the relative normalization of the fiare 
component in each annulus to scale with the area of the 
extraction region, allowing only its total normalization 
to be fit. Since the fraction of emission arising from flar- 
ing differed in each annulus, this procedure enabled us to 
disentangle its contribution more effectively than for our 
large-aperture analysis (reflected in better constraints on 
our parameters). Analogously to our large-aperture fit- 
ting, we initially allowed Zpc to be fitted freely and tied 
the remaining elements to Fe in their appropriate solar 
ratios. To begin, we tested the hypothesis of no temper- 
ature gradient (i.e. tying the gas temperature between 
each annulus). This isothermal model gave Zfo> 2.8, 
but was statistically unacceptable (x^/dof=260. 6/179). 
Freeing individual abundances did not enable an accept- 
able fit to be found, allowing us to reject the hypothesis 
of isothermality. 

We obtained a statistically significant improvement in 
the fit (at >99.99% likelihood, on the basis of an f-tcst) 
if we allowed the temperature in each bin to be deter- 
mined independently, although the fit was still unaccept- 
able (x^/dof=217. 6/175). Examination of the spectra 
revealed significant residuals at low energy, especially in 
the innermost annulus, analogous to those seen in the 
large aperture analysis. We therefore experimented with 
allowing the abundances of each metal in turn to be de- 
termined independently from Zpc, although significant 
improvements in the fit were only seen when O was freed 
(which fell to < 0.15). Adopting this model, we obtained 
a good fit to the data (x^/dof=189.4/174). We exper- 
imented with adding extra hot gas components to the 
fit to determine whether they led to significant improve- 
ments. However, we found no evidence of a significant 
change in either the abundances or the fit statistic. The 
best-fitting value for Zpo (shown in Table 0) was consis- 
tent with solar values (the lower limit at 99% significance 
on Zpe was 0.53). Furthermore, the gas was clearly defi- 
cient in O as compared to solar values; the 99% upper- 
limit on the Zo/Zpc ratio was 0.24. 

The temperature profile was relatively isothermal, ex- 
cluding a small peak in the centre (Fig g), although it 
was rather jagged in the second and third annuli. This 



is probably not a real effect; in fact kT in the outer two 
annuli is expected to be systematically high due to the 
truncation of the low-energy pass-band (Sect |4.5[ ) . Fur- 
thermore, due to the low temperature of the hot gas, the 
important O lines are blended to some extent with the Fe 
L-shell "hump" . Since the shape of the "hump" strongly 
constrains the temperature, this produces some parame- 
ter degeneracy between Zq and kT. As there were possi- 
ble systematic errors in our determination of Zq (Sect 0), 
we investigated the impact of such uncertainties by mea- 
suring the temperature profile with all elements tied to 
Fe. We found that the fit in this case was somewhat 
smoother, and the temperatures systematically higher 
(Fig. 0), although there was still evidence of a central 
temperature peak. 

As for the single-temperature, large-aperture analysis, 
we next systematically freed individually the remaining 
abundances and determined appropriate confidence re- 
gions for the parameters. Although we did not find any 
evidence of an improvement in the fit statistic by this 
method, this simply reffects best-fitting abundance ra- 
tios with respect to Fe close to solar (at least for those 
elements which could be constrained). We found that 
we were able to constrain Si, Ne and Mg in addition to 
O and Fe. In contrast to our large-aperture analysis, 
we were able to obtain reliable confidence intervals when 
all of these metals were simultaneously freed. The best- 
fit abundances are shown in Table 0. It is interesting 
to note that the abundance ratios were somewhat bet- 
ter constrained than the absolute abundance of Fe. This 
is easily understood since the primary source of uncer- 
tainty on the abundances is our constraint on the con- 
tinuum level. If the continuum is over-estimated, the 
equivalent widths of all the lines will be systematically 
under-estimated, and vice versa, producing a strong cor- 
relation between the abundances (manifesting itself in 
long, thin confidence regions in parameter space). We 
used a slightly modified implementation of the standard 
Xspec vapec model, in which the absolute abundance of 
Fe was directly determined, but for the other elements 
the fit parameters were the abundance ratios (in solar 
units) with respect to Zpo- This enabled us to measure 
directly the errors on the abundance ratios in Xspec with- 
out having explicitly to take account of thi s correlation. 
In previous studies (e.g. Buote et al. 2003b ), we adopted 
Monte-Carlo simulations to achieve this. 

For our 3D analysis, the infiated uncertainties in the 
profile introduced by the deprojection method were pro- 
hibitive if the metal abundances were to be determined 
by the fitting. We therefore fixed the metal abundances 
to the best-fit values obtained from 2D fitting and per- 
for mcd t h e depr ojection using the algorithm outlined 
in Buote ( 2000a ). To investigate the sensitivity of the 
temperature profile to our abundance determination, we 
additionally considered the case when Zo/Zfo= 1- We 
found the results in both cases were very similar to equiv- 
alent results found for our 2D analysis (Fig ^) . The tem- 
peratures tended to be a little higher in each annulus 
using 3D fitting, although this was most pronounced for 
the Zo/Zfc=1 case, which cannot acceptably fit the spec- 
tra all annuli for either 2D or 3D fitting. In contrast, the 
deprojected model with free Zq was able to give good fits 
in all annuli. The spectrum of the innermost bin and the 
best-fitting 3D model is shown in Fig |[ Since the tem- 
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Fig. 3. — Left panel: the temperature profile determined by fitting the projected spectra. We show results determined both with all 
elements tied to Fe (dotted-lines) and with Zq free (solid line). In addition, we show a best-fit power law model to the profile. Right panel: 
the same results, except fitting the deprojected spectra. For the deprojected analysis, it was necessary to fix Zpo and Zq values at their 
best-fit values. 



perature profile was relatively flat (excepting the inner- 
most bin), projection effects were not especially impor- 
tant. In the innermost bin, where there is a temperatm'e 
gradient, due to the relatively small spectral extraction 
region projection effects were also minor (Fig 0). 

4. SYSTEMATIC ERRORS 

In order to estimate the magnitude of any system- 
atic errors on our data-analyis, we considered explicitly 
a number of effects, which are discussed in detail here. 
Those readers uninterested in the technical details can 
proceed directly to Sect p. A summary of the magnitudes 
of each systematic effect is given in Table g. The dis- 
cussion herein focuses primarily on our simultaneous 2D 
spatially-resolved spectroscopy results, since they gave 
us the best constraints. 

4.1. Calibration issues 

A number of aspects of the satellite calibration may 
have possible systematic effects on our best-fit results. 
Using the standard CIA data-processing software it was 
not possible to correct S3 chip data for "charge transfer 
inefficiency" (CTI), which degrades the spectral resolu- 
tion of the data. We were, however, able to correct for 
a systematic drift in the instrumental gain by applying 
the apply_gain* task to our data An alternative proce- 
dure which allow s compensation f or th e effects of CTI 
was proposed by Townsley et al. ( ^002 )^, although the 
calibration for this method is far less up-to-date. Us- 
ing the CTI-corrected data (without the gain-correction), 
we accumulated spectra in a set of concentric annuli to 
match our spatially resolved spectroscopy. We obtained 
a good fit (x^/dof=184. 1/173) with the same spectral 
model, but with systematically lower overall Fe abun- 
dance (ZFe=0.79lo^26)- C)f the abundance ratios, only 

3-(-0.23^ 
-0.14i- 



Zo/Zpo was significantly altered (becoming 0.33_ 

We obtained a very similar temperature profile to that 

found for our standard analysis, although the central 



peak was less pronounced (kT=0.56to;o8keV). It is dif- 
ficult to disentangle to what extent the change in the 
best-fit abundances arose from the effects of the CTI 
correction, and to what extent they were simply a prod- 
uct of the less up-to-date calibration of this technique. 
Nonetheless, these values give an estimate of the calibra- 
tion uncertainty due to the failure to correct for CTI. 

It is also interesting to compare the results obtained 
for different releases of the Chandra Caldh; processing 
the data using version 2.25, for which the S3 spectral 
responses are known to have been subject to a small sys- 
tematic error led to systematically higher values of Zpe 
and Zo/Zpe (AZFe=i0.2 and AZo/Zpe^iO.lS). We can 
also use this discrepancy as a rough estimate of the mag- 
nitude of any possible calibration uncertainties remain- 
ing. 

4.2. The background templates 

A potentially serious source of error in the determina- 
tion of metal abundances is an uncertainty in the back- 
ground. Since the soft X-ray background is known to ex- 
hibit some dependence on the pointing, there were possi- 
ble systematic errors in our background estimation at low 
energies (the background template was found to match 
our observation well at high energies). In order to de- 
termine the extent to which background errors may be 
affecting our results, we simply experimented with alter- 
ing the background normalization by ±15%, comparable 
to the 10% variation between background fields found 
by M. Markevich^ . The best-fit metal abundances and 
temperature profiles showed no significant changes due 
to this procedure (Table 3), and the confidence intervals 
were not appreciably affected. Our results, therefore, 
were relatively insensitive to errors in the background 
estimation. 

4.3. The plasma codes 

In order to investigate the extent to which the choice 
of plasma code model can influence our results, we ex- 



* http://cxc.harvard.edu/cont-soft/software/corr_tgain. 1.0. html 

^ http://cxc. harvard. edu/cont-soft/software/ACISCtiCorrector.l. 37. htmi" see http://cxc.harvard.edu/contrib/maxim/acisbg/COOKBOOK 
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TABLE 2 

Abundance measurements and error budget 



Par. 



value AStat. Acalib Abkd Acode Aflare Abandv 



Asources AA^i 



H 



Zpe 


1.10 


Zo/Zpc 

Znc/Zpc 
ZMg/Zpe 
Zsi/Zpc 


0.05 
0.99 
1.07 
0.94 



+ 1.7 
-0.4 
+0.10 
-0.05 

±0.20 
±0.27 
±0.53 



±0.3 


±0.02 


-0.16 


±0.02 


±0.33 


-0.10 


±0.01 


+0.3 


±0.01 


-0.02 


±0.01 


+0.11 


±0.01 


±0.01 


±0.1 


±0.03 


-0.03 


±0.02 


±0.03 


±0.01 


±0.02 


±0.08 


±0.02 


+0.33 


±0.06 


±0.04 


+0.03 


±0.02 


±0.15 


±0.05 


+0.18 


±0.03 


±0.06 


+0.05 


±0.02 



Note. — Emission-weighted abundance measurements for the spatially-resolved spectroscopy, 
fitting a single-temperature model for the hot gas as a function of radius. In addition to the 
statistical (Stat.) errors, we present an estimate of the possible magnitude of uncertainties on 
the abundances due to calibration (calib), background (bkd), plasma code (code), flare model 
(flare), bandwidth (bandw.), unresolved source spectrum (sources) and Nn. It is difficult to 
assess the impact of multiple sources of systematic error simultaneously affecting our results 
(and it certainly is not correct simply to combine the errors in quadrature), so it is probably 
most correct to assume that the source of the largest systematic error dominates the systematic 
uncertainties. The magnitudes of the systematic errors listed are the changes in best-fit values, 
which were often more sensitive to such errors than the confidence regions (see text). 



perimented with replacing the APEC hot gas mode l with 



a MEKAL ([Mewe et al.| |1985| ; [Licdahl et al.| |l995D com 



ponent, since the atomic physics, particularly of the Fe 
L-shell line, is somewhat different. In general, both mod- 
els were able to fit the data, although MEKAL gave sys- 
tematically poorer x^ values (for example, the best-fit 
X^/dof for our best simultaneous fit to the spatially re- 
solved spectra was 202.2/174, as compared to 189.4/174 
for the APEC model). 

We consistently obtained Fe abundances which were 
~10-20% lower with the MEKAL code (and a sim- 
ilar effect with the Zq/Zfc value, although this was 
poorly-defined due to the large error-bars). The best-fit 
MEKAL abundances for this fit were Zfo=0.94+q42 and 
Zo/Zfo^O.OIIqJp Freeing up other abundances, we 
found good agreement with the APEC results for Ne 
and Si. However, ZMg/Zpc was systematically higher 
(1.43 ± 0.40), although marginal ly consistent within 



Buote et al.| ( ^003b| ) 



error-bars with the APEC result 

noted a systematic discrepancy of ^10-20% between the 
shapes of the APEC and MEKAL spectral codes, in the 
vicinity of the Fe L-shell lines. This effect resulted in 
^apcc ^ Zp^""^*^' for single-temperature fits, exactly as ob- 
served for NGC 1332. Since the Mg lines are close to the 
strong Fe L-shell lines (although, at the temperature of 
the system, they do not overlap) , there is probably some 
interplay between them in the fit. 

4.4. The flaring background model 

In order to obtain formally acceptable spectral-fitting 
results, it was necessary to include a term to account for 
the residual flaring, even though we filtered the data care- 
fully to remove the strongest flares. The flare was most 
evident above ~2 keV and had uniform surface bright- 
ness across the interesting regions of the S3 chip. Since 
the spectral shape of the fiare was not extremely well- 
constrained (due to the quality of the data) and since 
there were some discrepancies in its spectrum when esti- 
mated by two different means, this component was natu- 
rally a potential source of systematic errors. In order to 
determine the magnitude of any possible bias introduced 
by incorrect fiare modelling, we investigated the effect of 
replacing our preferred model with the alternative one 



based on our SI chip flare estimate. This gave a slightly 
worse fit (x^/rfo/=193. 7/174) but there were no signif- 
icant changes in the best-fit abundances, although the 
error-bar on Zpc was slightly widened so that the 90% 
upper-limit became 3.5. The best-fitting temperatures 
in each annulus were also not significantly changed by 
this procedure. 

As a further test, we fitted the data by omitting com- 
pletely any term to account for flaring. In the outer- 
most annuli (where the contribution due to flaring is 
most significant due to the larger extraction area), the 
model significantly underestimated the data. Since the 
spectral shape of the flare is signiflcantly harder than 
that of the unresolved point-source component, the un- 
resolved emission term could not move to compensate 
entirely for the flaring contribution, so that the overall fit 
was formally unacceptable {x^ /dof— 227/175). Nonethe- 
less, neither the best-fitting abundance values nor the 
confidence ranges associated with them were dramat- 
ically altered from our best-fit values (AZpc— ±0.3, 
AZo/Zpo— ±0.01). Examination of the residuals to the 
fit indicated that the model had simply fitted at low en- 
ergies at the expense of the data > 2 keV. The tempera- 
tures in all annuli were consistent, within error-bars, with 
the best-fitting 2D solution. However, the temperature 
was systematically lowered (by ~0.05 keV and 0.1 keV) 
in the outermost two annuli, respectively. Nonetheless, 
our abundance measurements appear relatively insensi- 
tive to the modelling of the flare. 

4.5. The bandwidth 

The ability of the spectral-fltting models to determine 
reliable metal abundances and temperatures is depen- 
dent to_somc_extcnt upon the bandwidth used in fltting 
(see Buote 2000a). In order to investigate the impact 
of narrowing our bandwidth, we first experimented by 
truncating the data above 2.0 keV. Such a procedure 
may be adopted to remove data heavily-biased by flar- 
ing, rather than the modelling we used. In adopting this 
procedure, we found that (omitting a flaring model com- 
ponent), we were able to obtain a good fit to the data 
(x^=158. 2/152). There were no significant changes in the 
best-fitting abundances or the temperature profile con- 
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sidering the size of the statistical error-bars. As might be 
expected, however, this procedure dramatically increased 
our error bars on the absolute metal abundances; the up- 
per limit on Zpo was no longer constrained, although the 
lower limit was not dramatically altered. The strong Fe 
L-shell lines in the 0.7-1.2 keV region play a vital role in 
determining both the temperature and the Fe abundance 
of the gas, and so it is unsurprising that we were still able 
to obtain a reliable gas temperature. However, since we 
did not have the higher-energy information, we were less 
able to constrain the continuum, hence the poorer con- 
straints on our parameters. 

In order to investigate the impact of our low-energy 
energy-cut, alternatively we truncated the data below 
0.7 keV and re-fitted the models. In this case, we found 
a downward shift in the best-fit iron abundance, but 
the predominant effect was to degrade significantly the 
error-bars on our results; we found that Zfc> 0.39 and 
Zo/ZFe< 0.61. The strongest relevant O lines are found 
at '^0.5-0.8 keV, so it is hardly surprising that removing 
the data at low energies reduced our ability to constrain 
the O abundance. The low-energy data were also critical 
in enabling us to constrain the continuum level, espe- 
cially as there were multiple components being fitted in 
each annulus. Failing to constrain the continuum trans- 
lates into poorly-determined metal abundances since it 
introduces serious uncertainties into the determination 
of line equivalent widths. Furthermore, the temperature 
profile was smoothed somewhat by this procedure, pre- 
dominantly through an increase in the best-fitting tem- 
perature in the third annulus (kT~0.38 keV) and a slight 
widening of the error-bars. Since we were forced to ex- 
clude data below ~0.7 keV in the outer two annuli during 
our analysis, this effect can explain the apparent "jagged- 
ness" of the temperature profile in these data-bins. The 
best-fit abundances will not have been as strongly af- 
fected by our truncating the bandwidth in these an- 
nuli (as is confirmed by the agreement between the two- 
temperature single aperture fit and the spatially-resolved 
analysis), since the abundances were tied between all the 
annuli. 

Since the background began to dominate the data at 
energies <0.5-0.7 keV, depending on the annuli, it was 
not possible to extend the bandwidth downwards in our 
annuli spectra to determine the impact of this proce- 
dure. However, for our large aperture extraction, it was 
possible to extend the bandwidth to 0.4 keV (although 
the calibration of the S3 chip is much more uncertain in 



this regime), as discussed in Sect 3.1. This increased the 



measured Zpc an d Zp/Zpc va l ues som ewhat (a similar ef- 
fect was seen by Buotc ct al. M03b), although the other 
abundance ratios were relatively unaffected. 

4.6. The unresolved source component 

We included a component to account for unresolved 
point-source emission in our data. The composite spec- 
trum of all the resolved point sources in the D25 re- 
gion of NGC 1332 could be fitted very well by a single 
bremsstrahlung model with temperature of 7.3 keV (see 
Paper I) . Although this is in excellent a greement with ob - 
servations of other early- type galaxies (Irwin et al. 2003), 
the spectra of fainter LMXB (which would be unresolved) 
tend to be somewhat harder than those of brig hter ob- 
jects (e.g. Church & Balucihska-Church 2001). This 
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Fig. 4. — The best-fitting spectrum to the innermost annulus 
(0-0.9'; 0-7.1 kpc) of the NGC 720 data. The data are shown, 
along with the best-fitting model. We also show the spectrum 
decomposed into each component; the hot-gas (APEC) model is 
shown as a dotted line, the unresolved point-source component 
(bremsstrahlung) is shown as a dot-dash line. The fiare component 
is shown with a dash-dot-dot line. 



discrepancy is likely to be exacerbated if high/ soft- 
state black-hole binaries are more common amongst the 
brightest source population of a galaxy. To determine the 
sensitivity of our abundance measurements to the precise 
shape of this component, we next allowed the tempera- 
ture of the bremsstrahlung term to fit freely. We found 
that its temperature tended to increase (although, due 
to the limited band-width of our spectra, it could not 
be constrained). The resulting temperature profile and 
abundances were in good agreement, within errors, with 
our best-fitting values. 

4.7. Column- density variation 

In our fitting, we fixed the neutral hydrogen column 
density to the canoni cal Galactic value es t imate d for the 
particular pointing (Dickey & Lockman 1990). In or- 
der to investigate any possible systematic effects arising 
from a mis-estimation of this value, we additionally fit- 
ted the data but allowed the iVn to fit freely. The best- 
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fitting column-density, Nii= (2.9;1;2 5) ^ 10 cm"', was 
in excellent agreement with our adopted canonical value 
(iVH= 2.2 X 10^° cm~^) and we did not find any signif- 
icant improvement in the fit statistic (Ax^=0.1) by al- 
lowing iVn to be free. Since our best-fit iVn was so close 
to the canonical value, we found that freeing Nn had no 
significant effect on the metal abundance or temperature 
profile determination. 

5. THE METALLICITY OF NGC 720 

NGC 720 is a very isolated, moderate-Lx elliptical 
galaxy with a slightly higher Lx/Lb than NGC 1332 and 
significant hot gas emission. Since it was observed by 
Chandra during AOl, it provides an excellent opportu- 
nity to investigate r nctal abunda n ces in another moder- 
ate Lx/Lb system. Buote ct al, ( 2002 ) reported Chan- 
dra analysis results for NGC 720, in which they were 
able to place tight constraints upon the shape of the 
dark matter halo. Unfortunately, there were strong peri- 
ods of flaring in the background during the observation, 
heavily contaminating the spectrum of the diffuse gas 
at energies > 2.0 keV. Since detailed spectral modelling 
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was not required in their analysis, these authors were 
able to mitigate against the effects of flaring by truncat- 
ing the spectrum above 2 keV. With such a truncated 
energy-band, it was, however, not possible to constrain 
the metal abundances (c.f. Sect p~^ ). We re-analyzed the 
Chandra data of NGC 720, adopting the techniques de- 
scribed in Sect. y| to estimate the spectrum of the flare 
from the S3 chip. The data were reprocessed identically 
to NGC 1332, and the periods of the strongest flaring 
were removed, giving a total exposure of 29.3 ks. The 
flare spectrum could be parameterized adequately us- 
ing the same model (broken power law plus broad gaus- 
sian line) as used for NGC 1332. Although the fit was 
not formally acceptable (x^/dof=234.4/196), the resid- 
uals predominantly appeared to be associated with low- 
amplitude, narrow features in the spectrum, so that the 
parameterization was a good approximation to the over- 
all flare shape. Since the flaring dominated the spectrum 
only over a small number of energy bins (at the highest 
energy), such features were relatively unimportant, and 
our failure to account for them simply degraded the x^ 
slightly. 

We extracted spectra from three annuli (with outer 
radii 0.9', 1.7' and 2.4') and centred at the X-ray cen- 
troid, computed in the same manner as for NGC 1332. 
The annuli were chosen to contain approximately the 
same number of photons, and extended out to 2.4'; we 
did not extend the annuli out further since the data 
were there very heavily contaminated by background flar- 
ing. We fitted the data simultaneously in the range 0.5- 
7.0 keV, with a model comprising a term for the hot gas 
(APEC), a term to account for unresolved point source 
emission (a 7.3 keV bremsstrahlung component), and the 
flare model. We found there was little difference in the 
quality of the fit if we allowed the temperatures to be 
determined separately or if they were tied between each 
annulus, and so we tied the temperatures to improve our 
abundance and temperature constraints. Similarly, we 
tied the metal abundances between each annulus so as 
to determine an emission-weighted average. We initially 
tied all the metal abundances to Zpc, in their appro- 
priate solar ratios, but we obtained a marginally unac- 
ceptable fit (x^/dof=400. 0/345). We next experimented 
with allowing each metal in turn to fit freely of Fe, and 
found a statistically significant improvement in the fit 
only if the O was fitted freely. We obtained a good fit to 
the data (x^/dof=368. 2/344) and the abundances were 
found to be ZFe=0.47J:o?^ and Zo/ZFc=0.07ig-^? and 
kT=0.55±0.02 keV 

Next, exactly as for the large aperture analysis of 
NGC 1332, we experimented with extending the band- 
width down to 0.4 keV, as the count-rate was still sig- 
nificantly in excess of the background in each annulus 
at these energies. Since NGC 720 was observed dur- 
ing AOl, the degradation of the quantum efficiency 
at low energies was a less significant effect and so 
the calibration of the S3 chip in this range should be 
more secure. As in NGC 1332, we found that our 
best-fit to the 0.5-7.0 keV band slightly overestimated 
the continuum between 0.4-0.5 keV (thereby leading 
to an under-estimate of the metal abundances). Fit- 
ting the model to the extended pass-band data, we ob- 
tained a good fit (x^/dof=383. 4/357) and somewhat 



higher metal abundances, but consistent with the error- 
bars of the restricted energy-range solution. We de- 
termined Zfo=0.71;°-^'j' and Zo/Zfo=0.23 ± 0.21, and 
kT=0.55±0.02. The lower-limit on Zfc, at the 99% con- 
fidence limit, was Zfc=0.47. We experimented with free- 
ing the remaining abundances, but were only able to ob- 
tain constraints on Ne and Mg (ZNe/ZFe=0.44 ± 0.43, 
ZMg/ZFo=l-26 ± 0.35). The strong flaring at high ener- 
gies largely swamped the data in the vicinity of the Si 
line, so that, in contrast to NGC 1332, we could not ob- 
tain Zsi/Zfo constraints. The best-fitting spectrum for 
the innermost annulus is shown in Fig. ^ These results 
were in general agreement with our results for NGC 1332 
(with the possible exception of the Zno/Zfc ratio, which 
was ~solar in NGC 1332, although the disagreement is 
marginal), and strongly excluded very sub-solar metal 
abundances except for O, for which Zo~0.2. 

6. DISCUSSION 

The Chandra observations of the normal, isolated, 
lenticular galaxy NGC 1332 enabled us to place con- 
straints on the metal abundances and the temperature 
profile of the hot gas. Revisiting the Chandra data of 
the isolated elliptical NGC 720 we were also able to place 
good constraints on temperature and the metal abun- 
dances. The temperature profile of NGC 1332 was rela- 
tively isothermal, excepting a slight temperature peak 
in its centre. For NGC 720, we found that the data 
similarly showed evidence of an '-^^isothermal profile. In 
neither case did we see convincing evidence of a central 
temperature dip, such as is seen in "cool core" clusters 
and groups. Both of these are normal, relatively iso- 
lated galaxies with moderate values of Lx/Lb, allowing 
us to investigate the enrichment of early-type galaxies in 
a regime between bright, gas- rich ellipticals (for which 
~solar abundances are typically being found; see intro- 
duction) and the very gas-poor galaxies in which very 
low abundances are still being reported. 

6.1. Near-solar Fe abundances 

In both of these galaxies, we were able to exclude the 
possiblity of the highly sub-solar Fe abundances often 
reported for early-type galaxies in the literature. These 
abundances were frequently subject to the "Fe bias" , on 
account of fitting a single tempe rature mode l to intrin- 
sically non- isothermal data (c.f. Buotc 1999). Multiple 
temperature components in a spectrum may arise from 
a strong temperature gradient over the spectral extrac- 
tion region or the mixing of photons due to projection 
effects (for 2D analysis). Alternatively, they may in- 
dicate a truly multiphase medium. We note that the 
hard component due to unresolved X-ray point sources 
in many early-type galaxies is too dissimilar spectrally to 
the hot gas to induce the "Fe bias" . The presence of more 
than one spectral component in the single aperture spec- 
trum of NGC 1332 can clearly be understood in terms of 
the slight temperature gradient (which is not so large as 
to require deprojected analysis) evident in our spatially- 
resolved results. This is supported by the good agree- 
ment between our single-temperature spatially-resolved 
and two-temperature single aperture abundances. The 
ability of of a two-temperature model to parameterize a 
temperature gradient in the observing aperture, and still 
yield reliable abundances, is an established effect (see 
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e.g. Buote et al. ^003b ). It is possible that the spectrum 
is somewhat more complicated than this picture. How- 
ever, the general agreement between the PLDEM and the 
two-temperature fits in the single aperture, in addition 
to the agreement between the spatially-resolved and sin- 
gle aperture fit results, suggest that the abundances are 
not especially sensitive to the precise parameterization, 



prov: ded it contains more than a single temperature com- 
po nent in this r egion. Similar conclusions were reached 



for the abundance ratios of several other elements. Even 
taking into account the likely systematic effects ob- 
served in our data, the low Zq/Zfo ratio we measure 
in both galaxies seems a robust result. Sub-solar val- 
ues of Zo/Zpc have been reported, using a variety of 
instruments, in a variety of other svstems su c h as the 
bright gro up NG C 5044 (iBuote et al.| | 2003b| : JTamura 
etal.||2003| ), M87 (IGastaldello fc Molcndi||200l ), the mas 



by Buote et al.| (p003§). 

For the spatially-resolved spectroscopy of NGC 1332 
and NGC 720, the lack of a strong temperature gradi- 
ent (which removed the need for deprojection) and the 
fact that the fits were not improved appreciably by using 
more than one hot gas component indicate that these 
data were not subject to the Fe bias. In fact, similar 
modelling of the ASCA data simply gave rise to very 
poor abund ance constraints (Zfo> 0.3 and Zfc> 0.2, re- 

to the so- 
Interest- 



sive elliptical NG C4636 (pCu et al.|p002| j, the starburst- 
ing galaxy M82 (Tsuru ct al 1997) and some clusters 



corrcctin; 



spectively: Buote fc Fabian 1998 

lar abundances of IGrevcsse fc Sauvall 199 



of galaxies ( Loewenstem 2004 , and references therein) 
We were also able to obtain constraints upon Mg, Ne 
and (for NGC 1332), Si. Intriguingly, we found that the 
Zo/ZMg ratios were also significantly sub-solar. These 
elements are primarily produced during Type II super- 
nova e, for which rnodels tend to predict Zo/ZMg ^0.9- 
1.9 ( [Gibson et al. 1997, convertin g to the solar abun- 
dances of [Grevesse fc Sauval 1998). Similar discrepan- 
cies have bee n found in several sy stems ( e.g. the centre 



ingly, their best-ht ASCA values (^1.2 and ^^0.6) were 
in excellent agreement with our Chandra results. Multi- 
temperature modelling of low-Lx systems with ASCA 
or Rosat freq uently gave similarly p oor abundance con- 
straints (e.g. Fabbiano et al. 1994), or were subject to 
parameter degeneracy in w hich erroneously low va lues of 
Zpe could be obtained (see Buote fc Fabian 1998). 

Our Chandra results, and similar res ults in the slightly 
lower -Ly /Lr radio galaxy NGC 1316 ( Kim fc Fabbiano 
^003 ) now reveal a self-consistent picture of ^solar abun- 



of NGC 5044 



dances from the centres of cool core clusters down to 
these moderate Lx galaxies. This is in contrast, how- 
ever, to the extremely sub-solar abundances which are 
still tending to be found in the lowest-Lx systems. If the 
low signal-to-noise of the data does not bias these low 
abundances, the lack of enrichment in the lower-mass sys- 
tems is a significant challenge to our understanding of the 
chemical enrichment process. Although the ^^solar abun- 
dances found in NGC 1332 and NGC 720 are far more 
consistent with more massive groups and clusters, we 
still did not obtain the several times solar values of Zpo 
typically pre dicted by enrichment models for individual 
galaxies (e.g. Giotti et al. 199l|; [Loewenstein fc Mathews 



1991| ; prighenti fc Mathcws||1999| J 



6.2. Abundance ratio implications for supernova 
enrichment 

A number of effects can enrich or dilute the hot gas in 
an early-type galaxy, such as sup ernovae, stellar mass- 
loss a nd gas inflow/ outflow (e.g. Matteucci fc Greggic 
1986 ). The a-elements are predominantly processed in 
supernovae and so the abundance ratios of these metals 
with respect to Fe provide a direct probe of the super- 
nova enr i chme nt history of the galaxy (e.g. Gastaldello fc 
Molendi 2002 ). Mass loss from intermediate-mass stars 
is expected to contribute significantly to the ISM of an 
individual galaxy. However, assuming no nucleosynthesis 
of the a-elements in the envelopes of intermediate-mass 
stars, the abundance ratios of the ejected gas simply re- 
flect the ISM make-up at the time the star formed. For 



enric pment dominated by very old stars, for example, 
stellar mass-loss will resemble type II supernovae, which 
were primarily responsible for the early enrichment. 
In addition to Fe, we were able to obtain constraints 



Buote et al 



2003b| ; M87: pastaldello 



Molendi|p002| ). In addition, in NGC 1332, the Zq/Zno 
ratio is also dramatically lower than predicted from Type 
II enrichment (~ 1.3- 1.9; Gibson et al. 1997), which ex- 
acerbates the discrepancy in the Zo/ZMg ratios. In 
NGC 720, the Zq/Znc ratio is closer to expectation, but 
the error-bars on the Ne abundance are large. 

Since low abundance ratios for O have been seen in a 
variety of galaxies with different instruments, it is un- 
likely to relate to calibration uncertainties or any bias 
due to an abundance gradient (which we would not be 
able to constrain). SNe metal yields are highly sensitive 
to theoretical modelling, so perhaps the simplest expla- 
nation for o ur results is a probl em with the SNe models 
considered. Loewenstein ( 2001 ) proposed that a signif- 
icant contribution of hypernovae to the enrichment of 
gas in rich clusters could explain similar low values of 
Zo/Zpc and Zo/Zsi, since the 0-burning region is signif- 
icantly enlarged in hy pernovae. However , the hypernova 
models considered by |Umeda et al. (2003) did not repro- 
duce very low values of both Zo/Zj^g and Zq/Znc, as 
seen in our data, although the authors acknowledge that 
yields of Ne and Mg are sensitive to the assumptions of 
the model. It may be of interest to note that recent 
Chandra observations of the old SNR N49B revealed SN 
ejecta significantly enriched in Mg, apparently with out 
accompanying Ne or O enrichment (Park et al. 2003). 

Notwithstanding these reservations, we were able to 
place some constraints upon the relative contribution of 
SN Type la and Type II to the ISM enrichment. In or- 
der to compare with pre vious results, we performed fit - 
ting analogous to that of Gastaldello fc Molendi (2003). 
Given the uncertainties in the hypernova yields, and to 
facilitate a meaningful comparison, we did not include 
any enrichment contribution from hypernovae. For a 
given fraction of enrichment arising from SNIa, f, theo- 
retical metal yields from SNIa and SNII ca n be compared 
with the data. We used the SNIa models of Nomoto et al 



( 1997b). and t h e wide lv-adopted SNII yields taken from 



Nomoto et al. ( 1997a , hereafter N97). In addit ion, we 



also c o nside red a range of SNII yields taken from Gibson 



ct al. (1997, hereafter G97), in order to examine the im- 
pact of the SNII model choice on our results. It is worth 
noting that these models all assume a Salpeter IMF. Al- 
though variation in the IMF is not expected to affect 
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the abu ndance ratios as siKn ificantly as the overaU abun- 
dances ( [Gibson et al. 1997), it remains a further source 
of systematic uncertainty. 

For the SNII yields of N97, we obtained / = 73 ± 5%, 
which was relatively insensitive to the SNIa model cho- 
sen (since, of the relevant elements, only Si is significantly 
produced in SNIa, and the Si error bar is large). For the 
range of SNII models considered in G97, we found sys- 
tematically lower enrichment fractions, reaching as low 
as / = 39 ± 10%. For NGC720, we found similar re- 
sults (/ = 85 ± 6% for the N97 model, and systemat- 
ically lower values if G97 yields were used, reaching as 
low as / = 62 ± 12%). Due to our large abundance ratio 
error bars, it was not possible to distinguish convinc- 
ingly between different SNIa and SNII models. One firm 
conclusion which we can draw, however, is that there is 
evidence of significant SNIa enrichment in both galaxies 
(f>0.4). Nevertheless, for comparison with other work, 
we adopt the N97 results as a standard, but caution that 
systematic discrepancies still exist between SNII models, 
which can significantly affect our results. It is interest- 
ing to compare our results with recent determinations of 
f ^80-90% for the centre of M 87, using the same yields 
(iGastaldello fc Molcndi| ^00^) . Similarly, in NGC 5044 
Buote et alj ( |2003bD fouSd f~70-80%, and f-80-90% 
was found in the centre of the bright nearby cluster 



A 1795 ( lEttorict al.||2002| ). The good agreement with our 
results further confirms consistency in the enrichment 
mechanism stretching from cluster-scales to moderate- 
Lx galaxies. 

Our fits to the data were, however, not formally ac- 
ceptable for any mixture of SNIa and SNII model yields 
(x^/dof> 4). Examination of the fit residuals revealed 
that, as expected, the low value of Zo/Zpe was prob- 
lematic, even if we allowed the value for NGC 1332 to 
increase by ~0.3, consistent with systematic uncertain- 
ties in its val ue. Uncertainties in the hypernova yields of 
Ne and Mg (Umeda et al, 2002| ) make it difhcult to in- 
vestigate similarly any putative hypernova contribution, 
which might account for the low Zq/Zfc ratio. If, in- 



stead, the low Zo/Zpc ratio reflects deficiencies in the 
SNII O yields, it may be sufficient simply to omit this 
data-point from the fitting. In this case we obtained good 
fits (albeit for only 2-3 data points), with a systemati- 
cally lower f for both galaxies. For NGC 1332, f=55±7%, 
using N97; using G97 it can be as low as 16 ± 13%. For 
NGC 720 we found f=68±8% for N97, and systematically 
lower values down to 25 ± 2 1 % for G97. In the centre of 
M 87 pastaldcllo fc Molcndi] ( |2002D also found a formally 
unacceptable fit; recomputing f based on their data, but 
omitting Zq/Zfc, we obt ained a good fi t with f=77ib5%. 
In the case of NGC 5044, [Buote et al'] ( ^003b[ ) computed 
f from the Si and S ratios, which were considered their 
most robust measurements, and so the low Zq/Zfc ratio 
did not similarly affect the results. Therefore, even when 
O was excluded from the fitting, we found good agree- 
ment between our results for NGC 720 (and, to a lesser 
extent, NGC 1332) and these more massive systems. In 
NGC 1332, the slightly lower value of f may refiect the 
increased significance of stellar-mass loss (from old stars) 
in a lower-mass system. 
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